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We present a derivation of the transformation between the Q and U Stokes parameters and the E and 
B scalar and pseudo-scalar fields. We emphasize the geometrical properties that such transformation 
must satisfy. We present the E and B decomposition of some simple maps and of a model for a 
supernova remnant. We discuss the relative amplitudes of the E and B components and argue that 
for generic random maps E and B should have roughly the same amplitudes. 



I. INTRODUCTION 

In the last few years there was a great surge of in- 
terest in the polarization anisotropies of the Cosmic Mi- 
crowave Background (CMB). The detection of CMB po- 
larization anisotropies has become a major goal in our 
field, prompting many groups to build experiments and 
to start thinking about future satellite missions dedicated 
to polarization (see for example |||-^). 

The pattern of polarization on the sky can be char- 
acterized in terms of a scalar (E) and a pseudo-scalar 
field {B) This decomposition is particularly useful 

because density fluctuations cannot produce B type po- 
larization . A B type pattern is a direct signature of 
the presence of a stochastic background of gravitational 
waves. Such detection would provide invaluable informa- 
tion about Inflation (for estimates of how constraints on 
parameters of the inflationary model would improve by 
measuring polarization see for example p|-pd|]). This is 
perhaps the most important source for the new interest 
in polarization. It has also been proposed that a detec- 
tion of B polarization could signal other types of "new 
physics" 11 . 

The mathematics of the E — B decomposition has been 
presented in several papers [p|J^Jl^-|l5[| . In this paper we 
will present a different derivation of the E — B transfor- 
mation that will highlight the ingredients that are needed 
to connect the spin two field of Q and U with the scalar 
and pseudo-scalar fields E and B. The aim is to gain 
intuition into the E — B decomposition, which is par- 
ticularly useful at this stage when new experiments are 
being designed. Intuition will help address issues such as 
the optimal shape of the sky patch needed to separate E 
from i?, or if both Q and U Stokes parameters need to 
be measured. 

We will also present the E and B decomposition of 
some simple polarized maps. Our aim is to understand 
if having a map with i? = 0, such as the one produced 
by density perturbation, is something generic or if one 
should always expect E k, B. 

The paper is organized as follow, in section §0 we 
present a derivation of the E — B decomposition, in §|l| 
we present the decomposition for some simple intensity 
and polarization maps. We comment about observational 
strategies and conclude in §[V. 



II. THE GEOMETRICAL PROPERTIES OF THE 
E~ B DECOMPOSITION 

The CMB anisotropy field is characterized by a 2 x 2 
intensity tensor lij. The intensity tensor is a func- 
tion of direction on the sky n and two directions per- 
pendicular to n that are used to define its components 
(§1,62). The Stokes parameters Q and U are defined as 
Q = (/ii — /22)/4 and U — /12/2, while the temperature 
anisotropy is given by T = {In -\- 122)/^ (the factor of 
4 relates fiuctuations in the intensity with those in the 
temperature, I cx T^). These three quantities fully de- 
scribed any state of linearly polarized light. While the 
temperature is invariant under a rotation in the plane 
perpendicular to direction n, Q and U transform under 
rotation by an angle '0 as: 



Q' ^Qcos2ip + Usm2il; 
U' ^ -Q sm2'ip + U cos 2-0 



(1) 



where e'^ — cos?/; §1 -I- simp 62 and §2 = — sin?/) §1 + 
cosV' §2. 

It is useful not to describe the polarization field in 
terms of Q and U but to do so in terms of two quan- 
tities scalar under rotation, usually called E and B . 
This E — B decomposition is a linear transformation of 
the Q — U field on the sky. The transformation in in- 
vertible. E and B differ in their behavior under a parity 
transformation, B changes sign while E does not. 

To make our derivation more transparent we will work 
in the fiat sky approximation, which is valid for small 
patches of sky. We do this only for the sake of clarity 
as all of our results can be trivially generalized to a full 
sky analysis. In the flat sky limit the directions (61,62) 
used to define the Stokes parameters at every point in 
the plane of the sky correspond to the coordinate axis, 
(61,62)= (x,y). 

A general linear transformation can be written as. 



B{e) = j (feKB{e,e)X{e). 



(2) 



where X is the 2 component vector X — {Q,U) and 



K 



iE,B: 



are the transformation kernels. 
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We want to derive the properties that the kernels must 
satisfy to make E and B transform correctly. We first 
consider two types of transformations, a translation and 
a rotation. Under a translation by a distance Oq the 
vectors on the plane of the sky transform as, 6' — O + Oq. 
Under a rotation of the coordinate system by and angle 
■0 they transform as, 6' — 'R.{ip)9, with R the standard 
rotation matrix. The explicit convention for the rotation 
is explained bellow equation (^. In both cases E and B 
should remain unchanged, in other words, 



E'{e') = E{e) (for translations 
B'{e') = B{e) k rotations). 



(3) 



Equation (ph implies that. 



j d2eKs(6/',e)X'(e) = j d^eKE{e,e)X{e) 
J (feKB{e',e)X'{e) = J (feKB{e,e)X{e) 



(4) 



Under a translation Q and U are remain unchanged un- 
changed, X'(0') = X{e). Equation (||) imphes that 
K(B,B)(0,e) = K(£;,B)(0-e). 

On the other hand, under a rotation Q and U are not 
scalars. They change as, 



x'{e') = iix{i')x{e) 



(5) 



with the matrix Rx(V') defined in equation (|^). Equa- 
tion (j^) implies that the E — B decomposition has to be 
non-local. E and B at point 6 cannot be constructed by 
combining Q and U at that same point because any such 
linear combination (if invertible) would not be scalar un- 
der rotations. 

The other type of transformation that needs to be con- 
sidered are reflections. After a reflection E remains un- 
changed and B changes sign. 



E'{e') = E{9) 

B'{e') = -B{e) (for parity). 



(6) 



Although equation (||) is valid for any reflection, to 
be concrete we consider a reflection about the y axis. 
The position vectors transform as 6' — {9'^,9'y) — 
{—Ox,Oy) and the Stokes parameters as Q'(9') = Q{9) 
and U'{9') — —U{9). The transformation laws for re- 
flections about other axis can be obtained by combining 
these transformation laws with the transformation prop- 
erties for rotations. 

Rather than trying to flnd directly the form of the ker- 
nels needed to satisfy all the above properties, for ped- 
agogical reasons we will use figures ^ and |2| to derive 
the kernels in a more intuitive way. We first consider 
the contribution to E{9) from a point e at a distance 
along y. We will assume that the contribution from this 
point in not zero. Note that with this particular config- 
uration the axis (§1,62) at point e are aligned with the 



vector e — 9. This is the reason we chose this set up. We 
will obtain the kernels for other configurations using the 
scalar nature of E under rotations. 

The contribution to E[9) from e, which we will call 
6E{9), is: 



5E{9) (xa%Q{e) + alU{e), 



(7) 



we 



^{E.B){d y) - {^Ie B)'"'{E,B) 



have introduced 
t'^7k r\)- E is invariant under 



where 
K 

refiections. In figure |1| we consider a reflection across the 
e ~ 9 line, the y axis. After this transformation 9' = 9 
and e' = e but the Stokes parameters change as, Q' = Q, 
U' = -U. This imphes that a.f = 0. 

To construct a quantity that is invariant under parity 
(E) , the contribution from a point separated by y can 
only involve Q. Our conclusion is a consequence of the 
particular geometrical set up but we will use the transfor- 
mation properties under rotations to derive the kernels 
for other directions. 



PARITY -> 



U = L-L 



Q'=i;-rrQ 

U' = I'-I'3=-U 



E' = E 



& 



FIG. 1. Contribution to E{0) from point e. The two 
panels are related by a parity transformation, a reflection 
across the y axis. The axis labelled by numbers are the ones 
used to define the Stokes parapeters. In the flat sky limit 
(61,62)= (x,y). 
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ROTATION -> 



U = L-L 




Q=i;-i', 
u'=r.-i'. 



e-0 



e'-B' 



E' = E 



6' 



FIG. 2. Contribution to E{9) from point e. The two panels 
are related by a rotation by an angle tp. 

We use figure || to understand what happens under ro- 
tations. When we rotate the coordinate system by an 
angle ipy the Stokes parameter are changed as described 
by equation (Q). The position angle of point e with re- 
spect to the X axis which we will call (f> changes from 
tt/2 to tt/2 + ip. We have to allow for the kernels to de- 
pend on this position angle otherwise E would never be a 
scalar under rotations given the transformation proper- 
ties of Q and U. Using the above argument about parity 
we concluded that, 



af (7r/2) = uj{9) 
af (7r/2) = 0, 



(8) 



where we have called uj{9) the value of the kernel at tt/2. 
We have explicitly included a dependence on the separa- 
tion because there is no reason why all the points along 
y should contribute equally. is a scalar so 



a 



,{tt/2 + ^)Q' 



6E'{6') 
-al{7:/2 + i,)U' 



5E{e) 



(9) 



This equation should be valid for arbitrary values Q and 
U . Combining equation (Q) and wc obtain, 



af(0) cos(2<^) 
Qf(0) =a;(0~)sin(2^). 



(10) 



In fact equation ([T^) is simple to interpret, only Qr = 
Ir — It, the difference between the radial and tangential 
intensities can be used to construct E and the weight 
can only depend on the distance 9. In other words 



when constructing E at we should use the radial and 
tangential unit vectors to define the Stokes parameters, 
{ei,e2) = {er,e^). In this frame only Qr contributes to E. 
We have proven that E{9) can we expressed as, 

E{e) = I (fe uj{e)Qr{e + e) 



(fe Lu{0) [Q{e + e)cos{24>) 



- U{e + 9) sin(2^)] 



(11) 



Any choice of weight w will produce a quantity that is 
scalar under rotation and invariant under parity. 

A similar argument can be used to show that the only 
way to construct B, a quantity invariant under rotations 
but that changes sign under reflections is, 

B{e) = I (fe uj{e)Ur 



(fie Lu{0) [Q{e + e)ain{24>) 



U{e + 9) cos(20)] 



(12) 



where now Ur is the U Stokes parameter defined with 
respect to the radial and tangential directions. In general 
LJ in equations ( pT| ) and (|l^) need not be the same. 

The usual definition of E and B corresponds to a par- 
ticular choice of the weight tj. 



uj{e) = -1/^2 ie > 0), 



(13) 



{lo{0) = but as will become apparent later this fact is 
not important for smooth fields). There are several rea- 
sons why this choice is made, some of which are easier to 
understand when working in Fourier space. As equations 
( pj] ) and (|l^) make clear, the E — B transformation is a 
convolution and thus becomes a multiplication in Fourier 
space. The choice of cj is such that the relation between 
E — B and Q — [/ is a simple rotation in Fourier space 
with no scale dependent factor. With uj given in equation 
( p^ ) the relation is. 



Q{i) 
U{1) 



[E{l)cos{2cpi) 
[£;(Osin(20O- 



- B{t)sin{2cpi)] 
B{l)cos{2<j,i)]. 



(14) 



Furthermore with this choice the ensemble average of 
P = + C/^ is the same as the ensemble average of 
E^ + B^. The sign convention on the other hand is chosen 
so that positive values of E generate a tangential pattern 
of polarization. The convention is rooted in the weak 
lensing literature which has identical mathematics and 
where the E field corresponds to the projected density k 
which produces tangential distortions when positive. 

For the purpose of finding a linear combination otQ—U 
that tests for the presence of gravitational waves or B 
type polarization any choice of uj is equally good. Other 
practical considerations such as the geometry of the ob- 
served patch of sky will probably be more important. In 
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weak lensing there is a long literature that deals with dif- 
ferent choices of w, to create for example measures of the 
enclosed mass that are more local than the 1/6'^ weight- 
ing. It is beyond the scope of this section to summarize 
that literature. 

We want the reader to take away three basic points 
from the above exercise: 

• The construction of E and B out of Q and U is by 
its very nature non-local. 

• To construct scalars under rotation at point 6 we 
need to average the Stokes parameters around cir- 
cles centered at 9 using the radial and tangential 
directions of this circle to define the Stokes param- 
eters {Qr,Ur)- The weight along the circle should 
be constant. 

• To construct E (a scalar) we need to average Qr 
and to construct B (a pseudo-scalar) we need to 
average Ur- 



III. EXAMPLES 

In this section we 

consider simple intensity-polarization maps to illustrate 
some of the properties of the E — B decomposition. As a 
byproduct we will understand better if the fact that den- 
sity perturbations produce only E type polarization is a 
unique prediction or if most sources of polarized emission 
have this characteristic. 



A. Simple Maps 

We start by considering a localized source of radiation 
of typical extent L and centered around the origin. We 
first intend to compute the E ~ B fields for points far 
away from this distribution. In this limit, equation jlT 
and ( ]l^ ) become, 



6*2 

cos(20) 
612 



(15) 



(feu{e). 



The choice of lo implies that the amplitude of E and B 
decay as 1/^2, What is more interesting is that whether 
E ov B are different from zero depends on the direction 
of observation. In fact there is no way to make B zero 
everywhere and keep E different from zero. E and B 
are only zero everywhere if the source is not polarized on 
average. 

In the case where the average polarization is along the 
X or y axis, E and B are, 



BiO) . a 



(16) 



where Q gives the average polarization of the source and 
ils is the solid angle it subtends. The direction depen- 
dence of E^B can be easily understood in terms of parity 
transformations. The mean polarization is invariant un- 
der reflections along directions where B — and changes 
sign along directions where B 0. 

As we discussed in the previous section, the E — B 
transformation has to be non-local. This implies that 
a localized source of emission will produce E — B even 
outside the region where Q and U are not zero. Our 
simple exercise has shown that if the source has a mean 
polarization, the typical size of the E and B components 
are the same outside the source. 

Let us now consider points inside the source. The first 
potential problem is that the weight function seems to 
diverge at zero distance. We consider a smooth Q — U 
field that can be expanded in Taylor series. We assume 
we are calculating E and B at point 9 and we expand Q 
and U around that point, 

Q{e + 6)^ q\q + e^Q^^iQ + ^MjQ^^jIq... 

u{e + 6)^ u\q + huAe + \o^3U,,,\q... (17) 

where we denote derivatives with respect to the different 
coordinates as ,i and summation of indecesis implied. To 
compute E and B we replace equation (|l^) into equations 
( prj ) and (12). The first terms to contribute are the sec- 
ond order ones because of the cos(2(/)) and sin(2(/)) factors 

in (0) and (|l|). The l/e^ in uj{e) is canceled by the f 
coming from the Taylor expansion. In fact there is an 
extra 6 from the d^O] there is no divergence at the ori- 
gin. We also conclude that E and B are most sensitive to 
the "quadrupole" pattern around 0\ the quadratic term 
in the Taylor expansion. 

As an example we can consider a filament as shown in 
figure ^ The emitting region has a length L along the y 
axis and a width I along the x axis, {L » I). We will 
assume that the Stokes parameters are constant along 
the filament and are zero outside. 

We now imagine doing the integrals in equations ( [ll| ) 
and (|l|) one circle at a time. As long as the radius of 
the circle is smaller than / the angular integrals cancel. 
It is an easy exercise to compute the values of E and B 
at the center of the filament, 



E = 


-c Q 


B = 


-c U 




fL/l 


c — 


V. 







dx 



(18) 



where x = 9/1. Equation ( p^ ) shows that E and B only 
receive contributions from rings larger than I and the 
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maximum contribution comes from 6 w y^3/2l. The con- 
tribution from far away rings are down by the 1/9 . 



to a great extent but some E and B are left. Which 
dominates at a particular depends again on the rel- 
ative orientation of the polarization and the separation 
vector between 6 and the center of the patch. Thus for 
a random pattern we expect similar levels of E and B. 




FIG. 3. Examples of polarization vectors inside filaments. 
The two circles indicate points that contribute with equal 
weight to E and B at the center of those circles. The con- 
tribution from points along the smaller circle cancel as one 
moves along the circle. The contribution from the second cir- 
cle is different from zero. In the case shown the contribution 
is mainly E. The first two filaments (labelled a), produce 
mainly E type polarization inside the filaments while the sec- 
ond two ( labelled b), produce mainly B type. 

Although we have presented a very simplified example 
we see two important features of the E — B transforma- 
tion. If the Q—U fields are constant over a scale Z, rings 
smaller than that do not contribute to E or B. Around 
a particular point, if the polarization vectors tend to be 
aligned or are perpendicular to the direction over which 
magnitude of the polarization is changing the pattern has 
a larger E than B. To have B the pattern has to form 
an angle of approximately 45° with that direction. In the 
case of a filament this is illustrated in figure |[ 

Finally we consider a case in which the polarization 
patter is very random, formed by regions of finite extent 
of typical size L inside which the polarization is constant. 
Different patches are independent. We sketch such a pat- 
tern in figure ^. We want to know how the typical values 
E and _B at a point inside a particular region compare. 
From our above examples we can conclude that the con- 
tributions to E and B coming from external patches are 
statistically the same, only depending for any particular 
external patch, on the relative orientation of the polar- 
ization in that patch and the separation vector 0. The 
contribution from points inside the same patch cancels 



FIG. 4. For illustration purposes we show a random pat- 
tern of polarization with a coherence length we call L. For the 
two points at the center of the filled circles, only the regions 
outside the circles contribute. External patches contribute on 
average the same to E and B. Only points inside the patch 
(but outside the circles) will contribute dominantly to E or to 
B. Which contribution is larger depends on the orientation 
of the polarization inside the patch and the position of the 
point where E and B are being calculated. 

Our final example argues that for polarization patterns 
that have a finite coherence length one should expect to 
have roughly the same E and B. This shows how remark- 
able it is that density perturbations do not produce any 
B modes. In order for B to be zero, the integral of Ur 
has to vanish identically (not just statistically) for every 
possible ring around any point, regardless of the radius 
of the ring. It is clear that this important symmetry will 
not hold for most random processes. 

B. Supernova remnant 

In this section we consider a more realistic model for 
polarized emission, a model for the emission of a super- 
nova remnant (SNR) It has been used successfully 
to model the emission from the the Galactic Spurs at 
radio frequencies around 1.4 GHz. 

The basic features of the model can be summarized as 
follows. Radiation is produced by synchrotron emission 
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from a shell. The thickness of the shell depends on po- 
sition and (to first order in the shell thickness) is given 
by: 



e{ip) — Ar/r — esin('0) 



(19) 



where r is the radius of the shell, is the angle relative 
to the direction of the initial magnetic field and e is the 
maximum thickness of the shell. 

The interstellar magnetic field had a strength Bq before 
the SN explosion. Later it is oriented along the surface 
of the shell and is amplified to a value, B = i?o/2e, a 
consequence of flux conservation. This model assumes 
that energy distribution of the particles responsible for 
the emission is of the form. 



N{E)dE = KE^^dE. 



(20) 



Moreover it assumes equipartition between the energy 
density in particles and magnetic field, 



K 



E-'<+^dE ' 



B^ 

Stt 



(21) 



Detailed derivations of these equations as well as pa- 
rameters that can fit different structures in our galaxy 
can be found in ||l^,|l^. It is not our objective to analyze 
what is expected from particular structures in our galaxy 
but rather to use this model to make a map of polariza- 
tion and compute its E' — _B decomposition to help build 
intuition. For this purpose the only two relevant param- 
eters are the angle (/?) between the plane of the sky and 
the unperturbed interstellar magnetic field (Bq) and the 
maximum width of the shell e. Without loss of generality 
we will assume that Bq lies in the x — z plane (the x — y 
plane is the plane of the sky). The projection on the sky 
of the unperturbed field is Bq cos(/3). 

We then compute the intensity and polarization ob- 
served along each line of sight be integrating the syn- 
chrotron emissivity along the line of sight. 



I = A / dxB 



?(7+l)/2 



Q = AU J dxB^2^^^^^ cos{2(f>) 

U = AIl J dxB^2^^'^'^ sin(2(/)), (22) 

where 11 is the degree of polarization of the synchotrom 
emission, B±^ is the component of the local magnetic field 
on the plane of the sky and (f) is the angle between B^ 
and the x axis and A is a normalization constant that de- 
pends on several parameters such as the number density 
of particles. 

Figure |^ shows the intensity and polarization map for 
the case (i = 45°. In figures ^ and ^ we show the cor- 
responding E — B maps. E and B extend outside the 
remnant. The intensity and polarization map has sym- 
metries of reflection across the x and y axis. The E map 



respects those symmetries while the B map does not. 
The B maps changes sign as one moves across the x and 
y axis. 




FIG. 5. Temperature and polarization map for a model 
of SNR with P = 45°. Rods indicate the magnitude of 




FIG. 6. E type polarization for a model of SNR with 

13 = 45° 
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FIG. 7. B type polarization for a model of SNR with 
/3 = 45° 

To understand the behavior of the Stokes parameters 
and the E and B fields we will look at one dimensional 
cuts along the x axis (perpendicular to the magnetic 
field). We show several examples in figure ^. Each of 
the columns corresponds to a cut at a different height 
along the remnant. 
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FIG. 8. One dimensional cuts of T, P = ^/Q^TTP, E and 
B across SNR models. Each of the three columns corresponds 
to a different height y — (—0.6,0,0.6) in units of the radius. 
In each panel we show the results for three different values of 
13 = (0,45°, 90°). The E curve for /? = 90° and y = has 
been divided by a factor of 10 so all curves in the panel are 
visible. 

From figure ^ we can conclude that: 

• Both E and B extend outside the remnant and de- 
cay as 1/9^. 

• The temperature and polarization patterns are in- 
variant under reflections across both the x and y 
axis. This implies that B is zero along both axis. 
E is invariant under reflections across the x and y 
axis but B changes sign. 

• Both E and B tend to peak at the edges of the 
SNR, where the T and P peak. We could think 
of the edge of the SNR as an example of the more 
idealized "filament" that we considered in the pre- 
vious section. For (3 = 0° and 45° we find that E 
is positive at the peaks, that means that the polar- 
ization direction is perpendicular to the direction of 
the "filament" . For f3 = 90° the situation is oppo- 
site and the polarization is parallel to the filament 
direction. 

9 E is larger than i? by a factor of a few in a "typ- 
ical" place inside the ring of maximum emission. 
In particular, at the edges of the SNR, this means 
that the polarization is typically like patters (a) of 
figure ^rather than patterns (b). 

• If /? = 90" B is zero everywhere because the the 
polarization pattern has reflection symmetry across 
any axis going through the origin. 

• If /3 = 90° the E component has a very large peak 
at the origin because the pattern is circular around 
that point. Note that there is no emission (T = 
P = 0) at the origin, where E has the maximum. 
This illustrates the fact that one cannot define a 
degree of polarization using E and B because their 
relation to Q and U is not local. 

It is important to realize that many of our conclusions 
follow from specific symmetry properties of the source. 
In reality this symmetries will be broken by real world 
complications such as inhomogeneities in the density or 
magnetic fields. As our examples in the previous section 
show, as the symmetries are lost the amplitudes of E and 
B become more similar. Moreover, because the E ~ B 
transformation is non-local, if only part of the SNR is in 
the observed field, the E and B decomposition will be 
different. 
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IV. DISCUSSION 

Detecting a B component in the CMB polarization 
field would be a great triumph for cosmology. As dis- 
cussed above the transformation between Q—U and E—B 
is necessarily non-local. Moreover, just from the geomet- 
rical requirements the only way to construct scalar and 
pseudo scalar quantities is to average Qr and Ur over 
circles. Thus the geometry of the patch of sky to be 
observed should be such as to allow for many different 
circles to be inscribed. 

The first generation of modern experiments, [jljjs) mea- 
sure the Stokes parameters in a ring on the sky. Thus 
only one circle can be constructed and so even if both Q 
and U are measured at every pixel along the ring, there 
is only one possible linear combination of the data that 
measures only E (the average of Qr along the ring) and 
one linear combination that measures only B (the aver- 
age of Ur along the ring) . All other combinations of the 
data receive contributions simultaneously from both E 
and B. In practice could not use these combinations 
because 1// noise make them unreliable. 

In general, one cannot construct circles around the 
points in the edge of the observed patch. Thus if the 
aim is to have the most possible linear combinations that 
are sensitive to either E or B but not to both, the best 
strategy is to make the observed patch of sky as round 
as possible. 

Another consideration is that for each point where E 
and B are calculated one needs to average either Q or U 
in the natural frame, the radial and tangential directions. 
This means that in order to be able to do it for as many 
circles as possibles one has to measure both Q and U in 
every pixel. 

The emphasis of our paper was to find linear combina- 
tions of Q and U that contain information about E or B 
alone. However it is perfectly possible to distinguish E 
and B type polarizations from the correlation functions 
of Q and U . Formulas that relate the power spectrum of 
E and B with the correlation functions of Q and U can 
be found for example in |^ . Distinguishing E and B this 
way does not rely on the shape of the observed region, 
as the correlation functions can be calculated just using 
pair of points. 

One should realize however that obtaining constraints 
on B from correlation functions comes at the price of 
larger error bars. It is clear that even though one can 
estimate the power spectrum this way, the errors in an 
estimate of B have contributions from both the power in 
E and B type polarization. Thus, because we expect the 
B signal to be smaller than the E one it is much better 
to directly find linear combination (rather than quadratic 
combinations) of the data that measure B. Suggestions 
of practical ways of separating E and B from correlation 
data have been recently presented in [18|. In reality a full 
analysis such as the one described in |15| will incorporate 
all the information available in a given experiment and 



should be preferred. 

We have also analyzed the E and B patterns expected 
for simple maps. We argued that in general one expects 
both E and B type polarization to have comparable am- 
plitudes although not necessarily equal. Whether E or B 
dominates at a particular place inside the source depends 
on the symmetries of the source. The E and B transfor- 
mation is not local so some E and B "leaks" outside the 
source, unless on average the source is unpolarized. The 
amplitude of E and B is the same outside the source 
on average but which dominates at a particular point 
depends on the relative orientation between the polar- 
ization direction and the separation vector. 
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